International Journal of Advanced 
Engineering Research and Science 


International Journal of Advanced Engineering Research and 
E Science (IJAERS) 
Peer-Reviewed Journal 
ISSN: 2349-6495(P) | 2456-1908(0) 
Vol-11, Issue-5; May, 2024 
Journal Home Page Available: https://ijaers.com/ eee 
Article DOI: https://dx.doi.org/10.22161/ijaers.115.1 7 


Alternative Algorithms for Time Series Water Wave 


Modeling 


Syawaluddin Hutahaean 


Ocean Engineering Program, Faculty of Civil and Environmental Engineering,-Bandung Institute of Technology (ITB), Bandung 40132, 


Indonesia 
syawalf1 @yahoo.co.id 


Received: 07 Mar 2024, 

Receive in revised form: 21 Apr 2024, 
Accepted: 03 May 2024, 

Available online: 11 May 2024 


©2024 The Author(s). Published by AI 
Publication. This is an open access article under 
the CC BY license 
(https://creativecommons.org/licenses/by/4.0/) 


Keywords— calculation 
conservation of energy. 


algorithm, 


I. INTRODUCTION 


Time series water wave modeling continues to evolve, 
particularly due to challenges in accurately modeling 
complex phenomena such as wave diffraction at breakwater 


Abstract— This research builds upon Hutahaean (2024a), focusing on 
time series water wave modeling, with substantial portions of this article 
derived from that prior study. Modifications were made to the calculation 
algorithm; specifically, vertical water particle velocity is now computed 
using the continuity equation, and water surface elevation is determined 
through kinematic free surface boundary conditions. While this revised 
algorithm enables the model to simulate wave breaking, it results in an 
excessively high breaker height and struggles to accurately simulate 
wave conditions post-breaking, particularly at large wave heights. To 
address these challenges, the kinetic energy conservation equation was 
integrated into the model. This addition facilitated the creation of a water 
surface elevation equation through the superposition of kinematic free 
surface boundary conditions and the kinetic energy conservation 
equation. This integration yields an equation that maintains a balance 
between changes in potential energy and changes in kinetic energy, 
enhancing the model’s capability to simulate both the shoaling-breaking 
process and subsequent wave conditions in shallow waters more 
effectively. 


Boundary Condition. This model is adept at representing 
shoaling but falls short in simulating wave breaking and 
post-breaking conditions in shallow waters. While some 
advancements have been made, as demonstrated by 
Hutahaean, S. Achiari H. (2017) and Hutahaean S. (2019), 


gaps. Traditional methods based on velocity potential 
theory often fall short in these areas. A more robust 
approach is the use of the Boussinesq equation, which has 
been significantly developed by researchers such as 
Boussinesq, J. (1871), Dingermans, M.W. (1997), Hamm, 
L., Madsen, P.A., Peregrine, D.H. (1993), Johnson, R.S. 
(1997), Kirby, J.T. (2003), and Peregrine, D.H. (1967, 
1972), among others. 


the ability to reliably simulate post-breaking conditions 
inconsistent, with some datasets accurately 
modeled and others not. 


remains 


Building on this knowledge, recent research by Hutahaean 
(2024b) has made significant strides. Utilizing velocity 
potential theory, this study identified that breaking 
characteristics are prominent in the Kinematic Free Surface 
Boundary Condition. Leveraging this insight, a new time 
series water wave model has been developed, incorporating 
the Kinematic Free Surface Boundary Condition for 
calculating water surface elevations. This model also 


The conventional time series water wave model primarily 
comprises two equations: the water surface elevation 
equation and the horizontal water particle velocity equation. 
The former is formulated by integrating the continuity 
equation—taking into account both water depth and vertical 
surface particle velocity—and applying the Free Surface 


utilizes the continuity equation to determine vertical water 
particle velocities. 
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The current model successfully simulates shoaling and 
breaking but tends to overestimate the height of breaking 
waves and fails to accurately model wave behavior as it 
progresses into very shallow waters. One significant 
limitation identified is within the Kinematic Free Surface 
Boundary Condition, which accounts for changes in 
potential energy due to variations in water surface elevation 
without a corresponding adjustment for kinetic energy. This 
oversight results in an energy imbalance, as there's no 
accounted source of energy to facilitate the observed 
changes in potential energy. 


To rectify this issue, the model has been refined by 
integrating the Kinematic Free Surface Boundary Condition 
with a conservation of kinetic energy equation. This 
integration ensures a balanced accounting of energy 
changes, which enhances the accuracy of the water surface 
elevation predictions. 


Additionally, modifications have been made to the Euler 
momentum conservation equation, as detailed in Hutahaean 
(2024a). The updated Euler Momentum Conservation 
Equation now strictly enforces that horizontal velocity 
changes occur only on the horizontal axis and vertical 
velocity changes only on the vertical axis. 


With these critical modifications, the enhanced time series 
model now more accurately simulates wave breaking, 
effectively extending its predictive capabilities to shallow 
waters near the coast. 


Il. WEIGHTED TAYLOR SERIES 


In this paper, we adopt the x — z coordinate system, 
where the x -axis represents the horizontal axis and the z - 
axis denotes the vertical axis. 


A weighted Taylor series is a truncated Taylor series limited 
to first-order terms, with the influence of higher-order terms 
encapsulated in coefficients known as weighting 
coefficients. 


For a function of two variables f = f(x ,t) 
f(x + 6x,t + 6t) = 


d d 
F(x, t) + Yet É + yrr L (D) 


Yt,2 and y, are weighting coefficients. 
For function f = f(x, z, t) the weighted Taylor series is 
f(x + 6x,z + 6z,t + ôt) 


= t) + ôt a 
= f (x,z,t) Yt3 dt 
d d 
+x + yz £ ..(2) 


Yt3> Yx andy, are weighting coeffcients. There is no 
difference between y, in f(x,t) and y, in f(x,z,t). The 
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basic values of these weighting coefficients are, y,. = 2, 
¥t3 = 3, Yx = Landy, = 1. 


The updated weighting coefficient values are presented in 
Table (1), wherein these coefficients are dependent on the 
optimization coefficient € . A larger € corresponds to a 
greater influence of higher-order Taylor series terms. The 
methodology for computing these weighting coefficients is 
detailed in Hutahaean (2023). However, the precision of the 
weighting coefficients listed in Table (1) surpasses that of 
Hutahaean (2023). 

In this particular time series model, the selection of 
weighting coefficients hinges on the wave amplitude. 
Specifically, a higher wave amplitude corresponds to a 
larger € , falling within the range of 0.15 < € < 0.35. 


Table (1) Weighting coefficients. 


E Yt,2 Yt,3 Yx Yz 
0.010 | 1.999797 | 3.004905 | 0.998792 | 1.011458 
0.012 | 1.999707 | 3.007157 | 0.998257 | 1.016713 
0.014 | 1.999600 | 3.009870 | 0.997625 | 1.023049 
0.016 | 1.999477 | 3.013062 | 0.996894 | 1.030512 
0.018 | 1.999336 | 3.016751 | 0.996064 | 1.039152 
0.020 | 1.999178 | 3.020955 | 0.995135 | 1.049022 
0.022 | 1.999002 | 3.025692 | 0.994107 | 1.060178 
0.024 | 1.998809 | 3.030982 | 0.992979 | 1.072679 
0.026 | 1.998599 | 3.036843 | 0.991751 | 1.086589 
0.028 | 1.998370 | 3.043296 | 0.990422 | 1.101976 
0.030 | 1.998124 | 3.050358 | 0.988994 | 1.118910 
0.032 | 1.997859 | 3.058049 | 0.987464 | 1.137468 
0.034 | 1.997576 | 3.066390 | 0.985834 | 1.157729 
0.036 | 1.997275 | 3.075398 | 0.984104 | 1.179778 
0.038 | 1.996954 | 3.085094 | 0.982272 | 1.203703 
0.040 | 1.996615 | 3.095495 | 0.980339 | 1.229598 


With the weighted Taylor series, the kinematic free surface 

boundary condition becomes, 
_ dn dn 

Wy = Yt,2 at + YxUn dx 

Wy = vertical surface water particle velocity 

U, = horizontal surface water particle velocity 


n(x, t), water surface elevation equation. 


Meanwhile, the total acceleration of horizontal and vertical 
water particles is, 

Du du yy duu du 

7 = ¥t3 apo ay YW, 
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Dw _ dw P dw P Yz dww 
dt 3 ae "ax 2 az 


II. DEPTH AVERAGE VELOCITY 
The model is formulated using depth average velocity as 
its variable, where according to Dean (1991), the horizontal 
depth average velocity is, 


t 
Daaa” dz 


U = horizontal depth average velocity 

u= the horizontal velocity integration coefficient 
D=h+yn 

D = total water depth 

h = water depth towards still water level (Fig. 1) 

7 = surface water level towards still water level (Fig.1) 


e| Aa water level 
x 
Zo 


U = u(x, Zo t) 


Fig (1). Depth average velocity concept 


Equation (4) was initially formulated to model very long 
waves, such as tidal waves, using a horizontal velocity 
integration coefficient 6,, = 1. In this study, the integration 
coefficient is calculated by defining the depth-averaged 
velocity as the velocity at a depth Z = Zọ below the still 
water level. Here, Zois negative, as shown in Figure 1. 
Consequently, Zo is defined in terms of the water depth A , 
specifically Zọ = —čh, where 0 < € < 1. 


The velocity potential equation solution of Laplace equation 
is 

o(x,z,t) = G(coskx + sin kx) coshk(h + z) sin ot 

G is wave constant 

k is wave number 


: 2T 
o is angular frequency, o = r 


T is wave period. 


a. Coefficients of Integration f„ and fw 
Horizontal water particle velocity from velocity potential is 


d 
u(x,z,t) = — a = —Gk(— sin kx + cos kx) cosh k(h 


+z)sinot 
Vertical water particle velocity is 
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d 
w(x,z,t) = = 
= —Gk(coskx + sin kx) sinh k(h 
+z)sinot 


Based on the definition of the horizontal depth average 


velocity, 

U = u(x, —éh,t) = —Gk(— sin kx + cos kx) cosh kh(1 
— €)sinot 

Thus, 


u__ coshk(h +z) 
U coshkh (1— €) 
Equation (4) is written into the equation for f,,, 


1 m 
Bu = UD oe dz 
1 n 
By = Dooshkh(1 a) kh 8) f cosh k(h + z) dz 
_ sinhk(h+n) 
Bu = kD cosh kh(1 — €) 
sinh kD 
Pu 


~ kD coshkh(1 — Æ 


In deep water, kD ~ kh = On, 0 is referred to as the deep 
water coefficient where tanh 0r = 1, and @ = 1.70 is used 
in this research 
Lo (5) 
uU Or cosh On(1-€) 

Although this equation is formulated in deep water, it also 
applies to shallow water, given the law of conservation of 
wave number (Hutahaean (2024b)) , where 
dk(h + z) 

dx mA 
Thereby k(h + z) is constant, unchanged against changes 
in water depth, as well as k(h + 7) is constant . 


The integration coefficient of the vertical depth average 
velocity Pw , is formulated in the same way, 


1 n 
Bw = owl dz 


1 7 
Bw = ee |S k(h + z) dz 
_ coshkD-—-1 
By = kD sinh kh(1 — Æ) 


B cosh 0m—1 
w ™~ @n sinh 62(1-€) 


b. Transformation coefficients 
In the derivation of equations for water surface elevation 
and water particle velocity, variables for surface water 


particle velocities, namely horizontal velocity u, and 


n 
vertical velocity wą, are introduced. To apply these 
equations more broadly, it is necessary to transform these 
surface velocities into depth-averaged velocities. 


Using the definition of the depth average velocity, 
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U, _ cosh k(h +n) _ cosh 61 
fun = y= cosh kh(1—€) cosh @n(1 — €) 


Or 


Uy = AynU 
Where the transformation coefficient a, is, 
cosh 61 
Qun = —— sean T) 


cosh 6m(1-€) 
Since U,U, and U,U,;U, is derived from u, Thus its 
distribution over space and time is the same as uy. 
Therefore, its transformation coefficient is the same as the 
transformation coefficient of uy. 
Unun = Aun UU 
UņUņnUn = uy UUU 


The vertical velocity transformation coefficient is 
_ Wn _ sinh 6m 
Qon = w Sahenae (8) 
A relationship is obtained, 
Wy = Awn W 
Wy Wn = Awp WW 


W pW Wy = Awn WWW 


nenen 
In this study 9 = 1.70, € = 0.32, are used, with 
Pu = 1.033, dun = 5.52, Ay, = 5.53 
The effect of 0 value is that the larger the @ value, the deeper 
the breaker depth and vice versa, the smaller the 0 value, 
the smaller the breaker depth, but has no effect on the 
breaker height value. 


IV. THE CONSERVATION EQUATIONS 
4.1.Weighted continuity equation. 
The continuity equation is formulated using the weighted 
Taylor series and by working on the conservation of mass 
principle,, 
k= + Yee = 0 tts (9) 
This equation is formulated under the condition that 
horizontal particle velocity only changed on the horizontal 
axis, as well as vertical particle velocity only changes on the 
vertical axis. 


4.2. Kinetic Energy Conservation Equation 

The kinetic energy conservation equation is formulated by 
reviewing the inflow and outflow of kinetic energy in a 
control volume. In the moving fluid mass contained kinetic 
energy Thus in the flow of the fluid mass there is also a flow 
of kinetic energy. The horizontal kinetic energy flow is 
u(pEķx) and the vertical kinetic energy flow is w(pE,,), 
where the kinetic energy on the horizontal axis is Ex, = 


w2 


2 
] while the kinetic energy on the vertical axis is E,, = 3g 


> 


p is the mass density of water. To formulate the kinetic 
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energy conservation equation, the control volume in Fig (2) 
is used. 


A 
ô 
Zh 
(7 
ô 
Zh 
Y f 
ôx ôx 
h2 h2 
Fig (2). Control volume. 
a. Space Averaging 


The velocity distribution on the sides of the control volume 
is not uniform, so an average velocity that represents the 
flow velocity on each side is required. 

Defined as 

fŒ, z, t) = uU(pEkx) 

h(x, z,t) = w(pExz) 

The horizontal kinetic energy input to the control volume 
through side 123 ata time t = t, is 


ôx ôZ 
nafga t) 
= fz) -n Z- y Z 


2 dx Z 2 az 


ô dx d 
hr = f(x- =z) = fz) -y 


_ ( ox ee t) 
f =f \x z7 z’ 


bx df 6z df 


[UGA Vic tia a, 


The kinetic energy input rate on the 123 side is defined as 


fiith th 
finput an 


Substitution of the equations of f,, f2 and fz obtained 


ôx 
input =h= f(x- .zt) 
= f(x,2,t) -y Z 


2 dx 


Using the same method, the following was obtained 


ôx ôxdf 
foutput =f; =f (x +z) = fz t) 1 Vx gx 
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ôz dh 


2 dz 
ôz dh 
h(x, z, t) + Yz az 


Oz 
Rinput =h,= h(x,z = 7’ t) = h(x, z,t) — Yz — 


6z 
houtput = hs = h(x,z + = t) = 


b. Time averaging 

The inflow and outflow of kinetic energy is observed at a 
very small time interval ôt. At that time interval, although 
very small, there is certainly a change in velocity. 
Therefore, an average velocity that represents the velocity 
at the time interval dt is required. Defined the average 


velocity at the time interval from t = — Z ke t = Z is, 

_ f(uze-¥) + f(x,z,t) + f(x,z,t +) 

f= 3 

Using the weighted Taylor series, the following was 
obtained 

fe = fz, t) 


Thus, for the inflow outflow process at an interval of ôt, the 
velocity at time t = t can be used. 


The kinetic energy inflow-outflow at time interval ôt is 


I0 = (firiput ~ foutput JZ ôt + (hinput ~ houtput )ôx ôt 
By working out the input and output equations on the sides 
of the control volume that has been formulated, the 
followings are obtained 


af ah 
10 = (- y,5% =) S260 $ (-r- ap) ost 


The Se causes a Mes in kinetic energy in the 
control volume, 

OE, = (Epy + Ekz) pbx5zZ 

The principle of conservation of energy, 


(OE xy + ÔEkz)pPôxőz = (- Vx OX = ôzôt 


+ (7, =) 8xôt 


The substitute the definition of f and h and work on the 
assumption of incompressible flow are 


(Egy + ÔEgz)pôxőz = (- ei 


+(- —y,p wpb bxdt 


i) 5250 


Both parts of the equation are divided by pôxôzôt, and 


worked out at very small dt close to zero, 
dE Kx dExz duEkx dwExz 


+ = -y, — py Seen. (10) 


at at X ax Z az 
This equation is the kinetic energy conservation equation. 


4.3. Euler’s momentum conservation equation in 
horizontal direction. 

By using the weighted Taylor series and by using the same 
fluid flow conditions as the continuity equation formulation 
where horizontal velocity only changes in the horizontal 
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axis and vertical velocity only changes in the vertical axis, 


Hutahaean (2024a) obtained Euler's momentum 
conservation equation in the horizontal direction is, 
a duu _ d TẸ dw 
Yrs at 2 dax ‘3ax) at 
= = (ww — ww) — ga er (11) 


V. VERTICAL WATER PARTICLE VELOCITY 
EQUATION. 
Vertical water particle velocity equation is formulated by 
integrating the continuity equation with respect to water 
depth. 


n 
ve Í ax 2 + YzWnņ — YzW-hn = 0 


This equation is written as an equation for w, where the 


bottom vertical water particle velocity is ignored, 


Integration of the right-hand segment is solved by Leibniz 
Integration (Protter, Murray, Morrey, & Charles, 1985) 


[ xe F =| F aB da 
ax T Ex Taz = lary 


TUN a dh 
Lo = d= | udz-u Uy a -Unay 


The integration of the first term of the right segment is 
solved by the depth average velocity concept and the bottom 
horizontal water particle velocity is ignored and substituted 
into the vertical velocity equation, 
Yx dUD dn 

oi =| tae hey) 
Transformation into depth average velocity equation, 

= —— (By, 2 = auu =) PEA (12) 
Where, 
D = h+nņ is the total water depth. The coefficient p4, 
defined in Equation (5). The coefficient ay, outlined in 


Equation (8), serves as the transformation coefficient from 
surface vertical water particle velocity to depth-averaged 
vertical water particle velocity. Similarly, @,,, specified in 
Equation (7) is the transformation coefficient from surface 
horizontal water particle velocity to depth-averaged 
horizontal water particle velocity. 


VI. WATER SURFACE ELEVATION EQUATION. 

Water surface elevation equation is formulated by using 
Kinematic Free Surface Boundary Condition and kinetic 
energy conservation equation. 


6.1.Kinematic Free Surface Boundary Condition. 


Weighted kinematic water particle velocity is 
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an an 
Wy = Yt2 dt + Yxun dx 
Written as water surface elevation equation, 
dn 1 dn 
a” mal - Yet 5) 
Added with depth average velocity variable, the water 
surface elevation equation becomes, 


an 1 dn 
are (ayn W — Var Gun U = 2) aan. (13) 


6.2. Integration of the conservation of energy equation. 
The energy conservation equation is multiplied by dz and 
integrated over the water depth, 


1 [ ye +f ara = 
2g \ Jn dt J, dt) 
1 d 
2 (=r, £1, ae — yw = w2,)) 


. 1 
Although there is an element of TA 


in both segments of the 


equation, it cannot be removed, this is to keep the unit of the 
equation the same as the unit of the Kinematic Free Surface 
Boundary Condition, which is m/sec. 


The integration is solved by Leibniz integration method and 
by working on the concept of depth average velocity and 
bottom water particle velocity is ignored, 

dn 1 ( dUU By dWW 


dt b dt 


Where, 


(e = Aun) oa + (Bw _ awn) n 
B,D 


The water surface change equation is the sum of the water 
surface change from the Kinematic Free Surface Boundary 


= 


Condition and the water surface elevation equation change 
from the energy conservation equation. 


dy 1 dn 
(1 tA = Pail Gun U =. N 
_ 2 (30u | Bwaww ausD_ gy y3 
A at E Bu dt taD 7 (Bu ax Gun U =) + 
Yzčwnņ 3 
ramy ) (14) 


Equation (14) is the final water surface elevation equation 
of the water surface elevation equation. In the right segment 
of the equation there is a change in kinetic energy which is 
the source of energy for the change in potential energy in 
the left segment, Thus this equation can be called the energy 
conservation equation, which is a balance equation between 
changes in potential energy and changes in kinetic energy. 
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VII. HORIZONTAL WATER PARTICLE 


VELOCITY EQUATION. 
Equation (11) is worked out on the surface at z = n, 
du, Yx dUyUy dn 
ns at 2 ax 9 ax 
By ate variable depth average velocity, 
Ve3@un se + SS = —g ...(15) 


VIII. NUMERICAL METHOD 
The space differential is solved by the Finite Difference 
Method while the time differential is solved by the 
predictor-corrector method. Details of the numerical 
methods as used by Hutahaean (2024a). 


IX. MODEL OUTCOME 
a. Wilson’s Criteria (1963) 


Wilson (1963) categorized wave profiles based on the ratio 
of wave crest elevation to wave height (Fig (3)). The wave 
profiles based on the comparison numbers are presented in 
Table (3). 


trough 


Fig (3). Water wave profile according to Wilson (1963). 


Table (3) Water wave profile criteria (Wilson (1963)) 


Wave type Nmax 
H 
Airy/sinusoidal waves < 0.505 
Stoke’s waves 0.505 — 0.635 
Cnoidal waves 0.635 —1 
Solitary waves =1 
a. Model execution at constant water depth.. 


In this section, the results of model execution at constant 
water depth (Fig (4)) are presented, with two water depths, 
h = 20.0 and h = 8.0 m. 
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n(0,t) = Asinot 
‘ 4 Still water level 


—> 


: h = constant 
Incoming wave 


sea bed 


Fig (4). Constant water depth 


As input, sinusoidal waves with the equation 7(0,t) = 
; ; 2 . ; 
Asinot, is used, where o = a T is the wave period, T = 


8.0 sec. While wave amplitude A is used from 0.3-1.3 m, in 
order to obtain the types of wave profiles at various wave 
amplitudes. 


a.l. Wave amplitude A = 0.30 m 

The model execution results at h = 20.0 m is shown in Fig 

(5), where 22% = as 
H 0.6 

sinusoidal profile criteria. The model execution results at 

Nmax = 0.37 = 

H 0.6 
0.62, which according to Wilson’s criteria belongs to 


Stoke’s profile. 


= 0.55, is quite close to Wilson's 


h = 8.0 m is presented in Fig (6), where 


On Stoke’s profile, there is a deflection at the transition 
from wave crest to wave trough at 7 ~ 0.0m. The 
deflection can be the difference between sinusoidal profile 
and Stoke's profile. In the sinusoidal profile Fig (5), the 
deflection is actually visible but still very weak, the 
deflection can be eliminated by reducing the wave 
amplitude, thereby at deep water depth h = 20.0 m, wave 
amplitude A=0.30 m 
sinusoidal profile to Stoke's profile. 


is the transition limit from 


n U&W 


70 80 90 100 110 
x (m) 
eta (m) ———U (m/sec) W (m/sec) 


Fig (5). Sinusoidal wave profile, A = 0.30 m, h = 20.0 m 
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a.2. Wave amplitude A = 0.60 m 


The model execution results at h = 20.0 m are presented in 
0.70 


Fig (7), where Tmax — —— = 0.583, which based on 


Wilson's criteria belongs to Stoke's profile. In addition, 
there is a deflection at the transition from wave crest to 
wave trough, this deflection is also found in the Stoke's 


profile in Fig. (6). Fig (8) presents the model execution 
results at water depth h = 8.0 m, where max = = = 


0.725, with the wave profile belonging to the cnoidal 
profile. Deflection changes to a small wave crest. 


0.4 
0.3 
0.2 
0.1 
0 
-0.1 
-0.2 
-0.3 
30 40 50 60 70 


n U&W 
Nmax = 0.37 M 


deflection 


x (m) 


eta (m) =——— U (m/sec) 


W (m/sec) 


Fig (6). Stoke ’s wave profile, A = 0.30 m, h = 8.0 m 


1 £ 
j=) 
0.5 (= 

z 
oğ X 
> ° : 
c / = 


`~ 
deflection 


60 70 80 90 100 110 
x (m) 


eta (m) =——— U (m/sec) 


W (m/sec) 


Fig (7). Stoke wave profile, A = 0.60 m, h = 20.0 m 


max = 0.87 m 


æ 
-0.5 deflection 
20 30 40 50 60 70 
x (m) 
eta (m) —— =U (m/sec) W (m/sec) 


Fig (8). Cnoidal wave profile, A = 0.60 m, 
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h=80m 


a.2. Wave amplitude A = 1.30 m 
The model execution results at h = 20.0 m are shown in Fig 


1.75 
(9), where 12% = == 
H 25 


criteria belongs to the Cnoidal profile. The wave height is 
reduced from H = 2.60 m to 2.50 m, that is assumed due 
to deflection in the form of small wave. In Fig (10), the 
model execution results at water depth h = 8.0m are 


= 0.7 which according to Wilson's 


presented, three wave crests are formed, the largest of which 
is the main wave. In this case, the main wave can be 
classified as a solitary wave with wave height H = 2.20 m, 
considering that the entire wave body is above the still water 
level. There is a reduction in wave height because some of 
the wave energy is used to form two small waves, one of 
which has a solitary wave 


profile. 
E 
1.5 D 
= 1 5 
03 
< 05 v 
e 0 Š 
-0.5 fun Fi i 
deflection 

-1 

20 40 60 80 100 
x (m) 

= cta (m) ———U (m/sec) == W (m/sec) 


Fig (9). Cnoidal wave profile, A = 1.30 m, 


h=20.0m 
£ 
3 F 
z ? i 
æ 1 R 
> 0 Š 
= Il 
-1 x 
-2 
10 30 50 70 
x (m) 
eta (m) ———U (m/sec) W (m/sec) 


Fig (10). Solitary wave profile, A = 1.30 m, h = 8.0 m 


The results from the model execution at a constant water 
depth revealed that sinusoidal wave profiles typically 
manifest in waves with small amplitudes. The formation of 
a wave profile is influenced by both the wave amplitude and 
the water depth. In deep water, a wave with a given 
amplitude might exhibit a sinusoidal profile, whereas the 
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same amplitude in shallow water could result in a Stoke's 
profile. Conversely, a Stoke's profile in deep water may 
transition into a cnoidal profile in shallow waters, and a 
cnoidal profile in deep water can evolve into a solitary 
profile when the water depth decreases. 


The distinction between a sinusoidal profile and a Stoke's 
profile is not only evident through the Wilson criteria but 
also by the characteristic deflection seen in the Stoke's 
profile, which appears as a flattened line at the transition 
from wave crest to wave trough. In contrast, in a cnoidal 
profile, this deflection manifests as small, gently sloping 
waves, and in a solitary profile, the deflection presents as 
two or more smaller waves. 


Additionally, the primary wave is often accompanied by a 
tail wave or multiple secondary waves. The phenomena of 
shoaling and breaking of these secondary waves, which will 
be further discussed in the following section, are critical 
aspects of coastal wave dynamics. 


b. Model execution at sloping bottom. 
In this section, the model is executed on a sloping sea bed 
(Fig (11)), with water depth 15.0-1.0 m, bottom slope 


14.0 


— = 0.093. 
150.0 


n(0,t) = Asinot 


© Still water level 
Incoming wave 
> 


150m 


Fig (11). Sloping sea bed 


In the simulation, a wave with a period of 8 seconds and an 
amplitude of 1.00 meter was input. Unlike previous 
depictions, the simulation results for this sloping bottom 
scenario are presented with the water depth h on the 
abscissa, allowing for a clear visualization of the 
relationship between water depth and wave characteristics. 
The waves transition from deep to shallow water, indicating 
that in the model, waves move from right to left. 


The simulation was conducted until the main wave reached 
a water depth of approximately ~ 1.3 m. The results are 
presented in Figure (12) and Figure (13). At this depth h = 
1.3 m, the height of the main wave stabilizes at 0.6 meters, 
and the wave profile transitions to a solitary profile, as 
shown in Figure (13). 
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24 breaking 1% breaking 


co EE 


2 


0 


n, U, W, CL 


1 6 11 
water depth h (m) 


eta (m) ——U (m/sec) 


——W (m/sec) 


CL (m) 
Fig (12). Wave condition after breaking. 


Notes: CL (Crest Line) is the line that connects the 
maximum water level elevation along the wave trajectory. 
In Figure 12, it is apparent that there are two breaking 
events; however, there are actually more than two, as will 
be demonstrated in the next section (Figure 14). The first 
breaking event, which occurs at greater water depths, 
involves the main wave, while the second event, occurring 
in shallower waters, involves the breaking of the tail wave 
or secondary wave. 


n, CL (m) 
o 
fo) wn me 


1 12 14 #16 18 2 
water depth h (m) 


eta (m) -——CL(m) 


Fig (13). Main wave profile at water depth ~ 1.3 m 


The wave condition after the first breaking in the surf zone 
is highly unstable, as observed in the unstable crest line 
(Fig. (13) and Fig. (14)). Wave instability in the surf zone 
is a well-documented phenomenon. Within this zone, the 
secondary wave experiences breaking, which is depicted in 
Fig.(14), where two breakings occur in two secondary 
waves. The wave profile is a solitary profile, as shown in 
Fig.(15). 
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n, U,W, CL 
WNPRFPORNWA 


3 3.5 4 4.5 5 
water depth h (m) 


eta (m) ——U (m/sec) 


—— \W (m/sec) 


CL (m) 


Fig (14). Breaking of the secondary wave. 


There are 3 breaking points with 3 breaking conditions, 
namely : 
a. Main wave breaking, H,, = 3.06m, 


hp = 5.17 m, “2 = 0.59 
hpr 


b. 1* secondary wave breaking 1 Hp, = 3.19 m, 


hp- = 3.87 m, = = 0.82 
hpr 


c.2™ secondary wave breaking Hp, = 2.3m, 

Apr = 3.45 m, = = 0.67 
Breaker height as suggested by Komar and Gaughan (1972) 
is: 
Hy, = 0.39 g'/5(TyH2)/s (16) 
Where Ty = 8.0 sec., Hy = 2.00 m, g = 9.81 m/sec? obtain 


Hp, = 2.46 m. The closest model result is the 2"4 secondary 
breaking wave Hp, = 2.2 m. 


Breaker depth index from Mc Cowan (1894), 
Hpr = 

ao 0.78 (7) 

The closest model result is the breaking of the 1st secondary 


H 
wave — = 0.82. 
hp 


P 


n (m) 


3.4 3.6 3.8 4 
water depth h (m) 


Fig (15). Breaking wave profile of the secondary wave. 
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Next, the breaking condition of the main wave is shown 
(Fig. (16)) with the wave profile being a solitary profile, 
Fig. (17). 


n, U, W, CL 
NPBPORNWA 


4 5 6 7 8 9 
water depth h (m) 


eta (m) ——U (m/sec) 


—— W (m/sec) CL (m) 


Fig (16). The breaking condition of the main wave. 


n, U, W 


5 5.1 5:2 5.3 
water depth h (m) 


eta (m) ———U (m/sec) W (m/sec) 
Fig (17). Wave profile of the main wave at the breaking 
point. 


The model results from the sloped bottom indicate that it 
can effectively simulate the shoaling-breaking process and 
subsequent movement in very shallow water depths. The 
simulation identifies three types of breaking waves: a 
primary breaking wave and two secondary breaking waves, 
with the wave profile adopting a solitary shape at the 
moment of breaking. If the simulation continues, there may 
be additional breakings, such as a fourth and fifth wave. 


Additionally, the model reveals a significant increase in 
vertical water particle velocity at the point of breaking. 
According to the continuity equation, a decrease in water 
depth leads to a greater disparity between input and output 
in the horizontal direction, which, in turn, escalates the 
vertical velocity. Thus, there is a likelihood that wave 
breaking may occur when the vertical velocity becomes 
excessively high. 
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X. CONCLUSION 


The study demonstrated that by incorporating Surface 
Kinematic Boundary Conditions as the water surface 
elevation equation and combining it with the kinetic energy 
conservation equation, a model was developed capable of 
simulating wave breaking until a reduction in wave height 
occurs in shallow water. This combination generates an 
equation that adheres to the energy conservation principle, 
establishing a balance between potential energy changes, 
indicated by variations in water surface elevation, and 
kinetic energy changes. 


The model, with modifications to its calculation algorithm, 
consistently generates four distinct wave profiles according 
to wave height and water depth: sinusoidal, Stoke's, cnoidal, 
and solitary profiles. Overall, the model effectively 
simulates water wave mechanics in both deep and shallow 
water environments. 


Future developments should focus on refining the 
estimation of breaker height and depth to align more closely 
with the results obtained from physical model studies 
conducted by previous researchers, as well as those derived 
from velocity potential theory. 
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